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In order for clarifying what are the essential thermal effects that govern the chiral phase transition 
at finite temperature, we investigate, in the real-time thermal QED, the consequences of the Hard- 
Thermal-Loop (HTL) resummed Dyson-Schwinger equation for the physical fermion mass function 
S_R. Since Eh is the mass function of an "unstable" quasi-particle in thermal field theories, it 
necessarily has non-trivial imaginary parts together with non-trivial wave function renormalization 
constants. The analyses so far completely dismissed this fact, despite it being one of the basic 
outcome of thermal field theory. The "approximation", which neglects the imaginary parts of the 
E_R, gives (non-trivial) constraint equations to be solved simultaneously, thus cannot be consistently 
carried out except in the trivial case, E_r = 0. In the present analysis we correctly respect this fact, 
and study, in the ladder approximation, the effect of HTL resummed gauge boson propagator. Our 
results with the use of numerical analysis, show the two facts; i) The chiral phase transition is of 
second order, since the fermion mass is dynamically generated at a critical value of the temperature 
Tc, or at the critical coupling constant ac, without any discontinuity, and ii) the critical temperature 
Tc at fixed value of a is significantly lower than the previous results, namely the restoration of 
chiral symmetry occurs at lower temperature than previously expected. The second fact shows 
the importance of correctly taking the essential thermal effect into the analysis of chiral phase 
transition, which are, in the previous analyses, neglected due to the inappropriate approximations. 
The procedure how to maximally respect the gauge invariance in the present approximation, is also 
discussed. 

PACS numbers: ll.lO.Wx, ll.15.Ex, ll.15.Tk, ll.30.Rd 



I. INTRODUCTION 



Although lots of efforts have been made to undestand the temperature- and/or density-dependent phase transition 
in thermal QCD/QED, we cannot have yet truly understood even the relation between the chiral transition and the 
confinement-deconfinement transition. Beginning of the relativistic heavy ion collision experiments at BNL-RHIC, 
which has been actively producing new results that may prove the generation of the quark-gluon plasma phase, has 
attracted an increasing interest in studying the physics in thermal QCD, thus has given us an encouraging time to 
proceed to further investigations of the mechanism of phase transition in hot and dense gauge theories, especially in 
QCD and QED. 

The Dyson-Schwinger (DS) equation is proven to be a powerful tool to investigate with the analytic procedure the 
phase structure of gauge theories, especially in the vacuum gauge theories [1,2,3]. However, we cannot say that, at 
finite temperature and/or density, the DS equation analyses of chiral and/or di-quark condensation have been carried 
out successfully. 

In the preceding DS equation analyses [4-8], the lessons from vacuum theories have been so simply applied to thermal 
theories without close examination. In most analyses the ladder approximation was used by simply neglecting all the 
hard-thermal-loop (HTL) effects [5,7], or only by taking improperly the HTL effects into the gauge boson propagator 
[6,8]. As a result they have missed the essential contribution of thermal gauge field theories, especially the important 
effect from the "dymanically screened" magnetic mode (having in general a momentum-dependent "mass" , though 
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being massless in the static limit). Many analyses, by fixing in the Landau gauge, ignored the fcrmion wave-function 
rcnormalization constants (WFRCs) by taking their naive tree values [5,6,7]. Furthermore, most analyses done in 
the real time formalism did not discuss the physical fermion mass function S^j itself of the retarded propagator 
[6,8], without any respect to the existence of imaginary parts, together with the inaccurate use of the instantaneous 
exchange (IE) approximation to the gauge boson propagation [6,8]. All such improper approximation methods may 
have caused the neglection of would-be-large contributions to the DS equation otherwise existed. 

Then we should seriously ask whether we could rely on the previous results of the DS equation analysis on the chiral 
phase transition as the real consequences of thermal gauge field theories. Considering the troubles in the previous 
analyses [4-8] mentioned above, we should make a re-analysis by studying the HTL resummed DS equation in the real 
time formalism, thus might giving a new understanding on the phase structure and the mechanism of phase transition 
in thermal gauge theories. 

Main interest of the present investigation lies in clarifying what are the essential temperature effects that govern 
the phase transition and also in finding how we can closely take these effects into the "kernel" of the DS equation. 
Essential procedures of our analysis can be summcriscd as follows; 

i) Firstly we use the real-time closed-time-path (RTCTP) formalism [9], and study the physical mass function Eij 
itself, not the En, of the retarded fermion propagator. 

ii) Secondly we correctly respect the fact that the is the mass function of an "unstable" quasi-particle in 
thermal field theories, namely that the E/j has non-trivial imaginary parts as well as non-trivial WFRCs. Neglection 
of imaginary parts and non-trivial WFRCs actually give (non-trivial) constraint equations to be solved simultaneously. 
This fact was totally dismissed in the preceding analyses. 

iii) Thirdly and most importantly, devoting our attension to closely estimating the dominant temperature-dependent 
contributions, we focus on studying the DS equation being exact up to HTL approximation: Both the gauge boson 
propagators and the vertex functions are determined within the HTL rcsummation [10,11,12.13], with which the gauge 
invariancc of the result at least in the perturbative analysis is guaranteed. With the HTL resummed vertex functions 
we can explicitly write down the HTL resummed DS equation [12]. 

iv) Finally in connection with ii) and iii) above, the gauge-paremeter dependent contribution must be carefully 
studied without fixing the gauge into some definite ones, such as the Landau gauge. 

The third point listed above is better to be taken step by step into the actual analysis of the DS equation. The 
advantage of the DS equation analysis lies in the possibility of systematic (step-by-step) improvement of the approx- 
imation to the integration kernel through the analytic investigation. 

Thus, in the present paper we present the result of our first step investigation in the strong coupling QED; focussing 
on what happens when we take into account exactly the HTL resummed gauge boson propagators, we investigate the 
consequences of the ladder (point vertex) DS equation with the use of the fully HTL resummed photon propagator. In 
this approximation, however, the gauge invariance of the results is spoiled. The procedure how to maximally respect 
the gauge invariance of the results in the ladder approximation, is discussed and the result is given in a separate 
paper [14]. Analysis in QCD and effects of fully including the HTL resummed vertices will also be presented in the 
separate paper [15]. 

This paper is organized as follows. In the next section II, the DS equation in QED for the fermion self-energy in 
the HTL approximation is derived, and the improved ladder (point vertex) DS equation with the HTL resummed 
photon propagator is given. The instantaneous exchange (IE) approximation for the photon propagation is disccused. 
In section III wc; give the results of the present investigation with the use of numerical analysis: the order of the 
chiral phase transition, the critical curve in the T-a plane (where T is the temperature of the environment and a 
the coupling constant), and values of the critical exponents. Conclusion of the present analysis and discussion on the 
results are given in the last section IV. 

II. THE DYSON-SCHWINGER EQUATION IN QED FOR THE FERMION SELF-ENERGY FUNCTION 

A. The DS equation in the Hard-Thermal-Loop approximation 

The DS equation for the physical, i.e., the retarded fermion self-energy function E^ in the HTL approximation can 
be obtained by applying the following approximation to the full DS equation; 

i) replace the full gauge boson propagator with the HTL resummed propagator, and 

ii) approximate the full vertex functions to the HTL resummed vertex functions. 

Then in the RTCTP formalism we get in QED the desired DS equation (in QCD, the running coupling should be 
used inside the loop-momentum integration) [12], 
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X [*^RAAi-P^K,P-K)SRA{K) *T-^iAA-K,P,K-P) *Grr.,^4P - K) 
+ *T^raa{~P^K,P-K)Srr{K) *T''aar{-K,P,K-P) *Gra^^,{P ~ K)] , (III) 

where *G'"^ is the HTL resummed gauge boson propagator, whose retarded R = RA component is given by 

* (T 1 1 i 

and the RR component by 

*G''^^{~K,K)^{l + 2nB{k,))[*G''^{K)- *G^^{K)], ^^(fco) = --^^-L__, (11.3) 

with *n|? and *n]^ being the HTL contributions to the retarded gauge boson self-energy of the transverse and 
longitudinal modes [16], respectively. A'"^, RP'^ and D^'^ are the projection tensors [17], 

RPK" RPR" 
AP" = gP" - RP" - DP", RP" = ^, DP" = ^, 

y ' R^ ' R^ ' 

where R = (fc, kok), k = -s/k^ and k = k/Zc denotes the unit three vector along k. 

S{—P, P) = S{P) denotes the full fermion propagator, whose retarded RA = R component is given by 

Sr{P) = jp^.^ ^ , (11.4) 

and Srr the RR component 

Srr{P) = (1 - 2nF{po))[SR{P) - Sa{P)], npipo) = T^FYXT' ("-5) 



The fermion self-energy function S/j in Eq.(II.4) can be tensor-decomposed as 

Er{P) = (1 - A{P))p,Y - S(P)7° + C(P) , (II.6) 

with A{P), B{P) and C{P) being the three independent scalar invariants to be determined. At zero temperature, the 
wave function renormalization constant A{P) coincides with B{P) and equals to unity in the Landau gauge, while at 
finite temperature it is not. G{P)/A{P) plays the role of mass function, in which we are interested. 
*r'^ is the HTL resummed 3-point fermion-gauge boson vertex function [12], 

*Kp, = 7^p, + SKp,^ (a,/3,7 = A,i?), (H.7) 

where ST'^^^ represents the HTL contribution to the 3-point fermion-gauge boson vertex function, and 7^^.^ the 
tree vertex with Jraa ~ ^aar. ~ 7^' otherwise zero. Appearence of the HTL resummed vertex functions together 
with the HTL resummed gauge boson propagators assures that the HTL approximation is consistently carried out 
in studying the HTL resummed DS equation [12], and guarantees the result being gauge invariant, at least, in the 
effective perturbation regime. 



B. The improved ladder DS equation with the HTL resummed gauge boson propagator 



Neglection of ST'^^^ in Eq. ( II. 7 ), simply brings us to the ladder (point-vertex) DS equation with the HTL resummed 



gauge boson propagator. It significantly simplifies the structure of the DS equation to be examined, thus reducing the 
technical difhculty in handling the DS equation itself. The price to pay is to lose the assurance of gauge invariance of 
the results. 

In this paper as already mentioned before, we investigate the consequences of the ladder DS equation with the 
HTL resummed gauge boson propagator, Eqs. ( p^I.2| ) and (II. 2). In this case three invariants A{P), B{P) and C{P) 
as functions of po and p satisfy the following coupled integral equations {PP — (pojP))i 
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-p'[l-AiP)] = 



-B{P) 



C{P) 



(2^)4 



{1 + 2nB(po - fco)}/m[ *Gg'(P - K)] 



{K„Pp + ifpPcr - Po{K„gpo + Kpg„o) - kQ{P„gpo + Ppgao) + pkzg^p 

A{K) 



-'^PofJaOgpo} 



[ko + B{K) + ie]2 - AiKYk"^ - C{KY 
ko + B{K) 



[fco + B{K) + ieY - A{K)^k^ - C(KY 



+ {1 - 2nF(fco)} X 



GP^{P - K)Im {K^Pp + KpP, - po{K,gpo + Kpg^o) - kgiP^gpo + PpQao) 



+pkzgap + 2po^o5<TOffpo} 



A{K) 



'{PagpO + PpQaO - 2pog^Qgpo} 



[ka + B{K) + ie]2 - ^(AO^P - C{KY 
fco + B{K) 



(2^)4 

{Kagpo + Kpg^o - 2kog^QgpQ} 
-{25po25<TO - 9ap} 
*GP^ {P - K)Im 
-^kog^ogpo} 



[fco + B{K) + ie]2 - A{K)^k^ - C{KY 

pa 
R 



{1 + 2nB(po - fco)}/m[ *Gg'(P - X)] x 



A{K) 



[fco + + ie]2 - A(A:)2fc2 - c{KY 

fco + ) 



[fco + B{K) + ie]2 - ^(i^)2fc2 _ C{KY 

MK) 

[fco + B{K) + ie]2 - A{K)^k^ - C{K)^ 
fco + BiK) 



+ {1 - 2nF(fco)} X 

{KagpO + -fS^pffo-O 



[fco + B{K) 



- AiKYk"^ -C{KY 



{2gpo2g„o - g^p} 



+ 2nB{p, - fco)}/m[ *G''^{P ~ K)] x 



C{K) 



[fco + + ie]2 - A(A')2fc2 - C(if)2 



+ {1 - 2nF(fco)} X 



*G''^{P - K)Im 



C{K) 



[fco + B{K) + ie]2 - A{KYB - C(A:)2 



where the cutoff scale A is introduced in order to regularize the integrals over fco and k. 



(II., 



(II.9) 



(11.10) 



C. The ladder DS equation in the improved instantaneous exchange (IE) apporoximation 



The DS equations, Eqs. ( |II.8| )-( II.10 ), are still quite tough to be attacked, forcing us further approximations for 
the analysis to be effectively carried out. However, the approximation made use of must be consistent with the HTL 
approximation, without missing the important thermal effects out of the kernel of the DS equation. 

Here it is worth noticing that the instantaneous exchange (IE) approximation to the gauge boson propagation, 
frequently used in the preceding analyses [6,7,8], is not compatible with the HTL approximation in a strict sense. In 
the exact IE-limit the HTL resummed transeverse mass function, *Iij^{K), vanishes and the transeverse (magnetic) 
mode becomes totally massless. Namely the IE approximation discards the important thermal effect coming from the 
Landau damping, thus dismissing the dynamical screening of the magnetic mode. This causes the famous quadratic 
divergence of the Rutherford scattering cross section. 



To see this point more clearly, let us take the IE-limit of the DS equation, Eqs. ([I.S)-( II.IO ), and neglect Im[A{P)\ 
and Jto[C(P)], then we obtain the following equation for Im[B{P)], 



Im[B{P)] 



47r » 



k^dk 



dz 



E WE"^ 



]2 ' 



,212 



(11.11) 
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TABLE I: Sample values of invariants, Re[A{P)], , Im[C(P)], at po = 0, p = O.IA for a = 4.0. 



T Re[A{P)] Im[A{P)] Re[B{P)] Im[B{P)] Re[C{P)] Im[C{P)] 



0.130 


1.262 


0.125 


-0.547 


0.303 


0.496 


-0.035 


0.135 


1.227 


0.067 


-0.551 


0.311 


0.466 


-0.057 


0.140 


1.248 


0.008 


0.496 


0.347 


-0.341 


-0.071 


0.145 


1.249 


0.016 


0.416 


0.342 


-0.157 


-0.042 


0.150 


1.217 


0.059 


0.378 


0.315 


0.000 


0.000 



showing Im[B{P)] to be quadratically divergent. In order to ignore the Im[B{P)], we must disregard the divergence 
and set oo = 0. This fact means that we cannot naively neglect the imaginary parts of the S^j without facing the 
inconsistent constraint equations. 

The reason why in the previous analyses this divergence did not appear, is that the imaginary part of was 
completely neglected there from the beginning, namely that the equation for Im[I]fi] was totally disregarded. As is 
evident from the above, however, the "approximation" made use of in the previous analyses [6,8] is not the consistent 
approximation. 

Taking the above into account, the approximation we further make use of is the improved IE approximation to the 
longitudinal gauge boson propagator, by keeping the exact HTL resummed transeverse propagator. In the IE-limit 
the HTL resummed longitudinal mass function, *Ilf^{K), has a definite thermal mass ~ i^T)'^, representing the 
Debye screening due to thermal fluctuation, thus even in the IE limit the longitudinal mode can take into acount the 
essential thermal effect. In the present analysis the gauge is fixed to the Landau gauge = 0). 

Here it is fair to note that in the point vertex ladder approximation, as already mentioned before , th e gauge 
invariance of the results is spoiled. To maximally respect the gauge invariance, we should solve Eqs. (|L8|)-( |II.10| ) 



with the constraint A{P) = 1, which guarantees Z2 = 1, being consistent with the Ward i dent it y Zi = Z2. This can 
be done by successively adjusting the gauge-parameter ^ in solving the above equations ([I.8)-( II.1C ). The result of 
this gauge-parameter dependent analysis will be published in a separate paper [14]. 



III. RESULTS OF THE NUMERICAL ANALYSIS 



A. Numerical solutions 



Now we should solve numerically the DS equations Eqs. ([I.8)-( [I.10D with the IE approximation to the longitudinal 



mode [18]. As we have already stressed, we respect the existence of the non-trivial imaginary parts, /r7i[yl], Im[B] and 
Im[C], together with the corresponding real parts, Re[A], Re[B] and Re[C]. We use the numerical method consisting 
of starting with suitable trial functions for the solution and iterating the calculational procedure until stable solutions 
are obtained. This method is simpler to handle and is useful so long as the convergence of the iteration is guaranteed. 

We choose several types of the trial functions: various choices of constants, independent of po and p. 

At each iteration, the three-fold integration is performed, namely over k, over /cq, which are cut off at the mass scale 
A, and over z = cos 9. The integration kernel of the present DS equation shows a little bit singular behavior, and the 
numerical integration of such a singular integrand needs careful integration prescription, which is properly managed. 
As a result, in each calculation we performed at most 1000 times iterations, and obtained fairly stable solutoins. 

Result of the present analysis shows that the wave function renormalization constants receive 10-20 percent correc- 
tions, see Figure 1 and Table I, indicating the necessity of gauge-parameter dependent analysis [14]. The generated 
size of the imaginary part is nearly the same as that of the real part, showing the existence of the non-trivial imaginary 
parts, which is as expected though completely neglected in the previous analyses. 



B. Chiral phase transitions 



Let us now show the behavior of the mass fimction M (P) = Re[C {P) / A{P)] at some fixed value of p as a function 
of the parameters a = e^/(47r) and T. All the data shown in this section are those calculated at po = 0, which 
is suitable for studying the static mass. The T-dependence of the mass function M{P) = Re[C {P) / A{P)] with 
p = O.IA is shown in Figure 2 for various fixed values of a, and in Figure 3 the a-dependence of the mass function 
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FIG. 1: Typical shape of the function Re[A{P)] at po = for a = 4.0. 



M{P) = Re{C {P) / A[Py\ with p = O.lA is shown for various fixed values of T. The errors resulting from the 
fluctuations are smaller than the size of the symbol used for each sample point in Figures 2 and 3. 
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FIG. 2: The T-dependence of the mass function M{P) = Re[C (P) / A{P)] with p = O.lA for various fixed values of coupling 
constant a. 



Prom these figures we can see the two facts; i) The chiral phase transition is of second order, since a fermion mass 
is generated at a critical value of the temperature T or at the critical coupling constant a without any discontinuity, 
and ii) the critical temperature Tc at fixed value of a is significantly lower than the previous results [6,7,8], namely 
the restoration of chiral symmetry occurs at lower temperature than previously expected. The second fact shows the 
importance of correctly taking the essential thermal effect into the analysis, which was disregarded in the previous 
analyses due to the inappropriate approximations. 
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3: The a-dependence of the mass function M{P) = Re[C{P)/A{P)] with p = O.IA for various fixed values of temperature 



C. Critical curve and critical exponents 



Analyzing those data, typically shown in Fignrcs 2 and 3 in the last subsection, we can draw the critical curve 
that shows the phase-boundary in the (T, a)-plane. The critical curve and the critical exponents can be determined 
systematically from the data with the use of the following method. 

Wc fit the mass function M{P) with n data points, Mj, aj and Tj, with i = 1,2, ....n near the critical point by 
assuming the functional form 



M(P) 



for fixed coupling constant a, and 



M{P) 



'{a 



(III.1) 



(III.2) 



for fixed temperature T. Here, Ct, Ca, T^, ac, v and r] are adjustable parameters, with and corresponding to 
the critical temperature and critical coupling constant, and v and t] designate the critical exponents. The values Ct, 
Ca, Tc, etc, and "q can be determined with the method of the standard least-squares fit. 

The critical coupling constant a,, as a function of T, and the critical temperature as a function of a, thus 
determined are shown in Figures 4 and 5. By putting them together we can determine the critical curve in the 
(T, Q!^^)-plane, Figure 6, which shows two characteristic behaviors: a) As T becomes smaller, the critical coupling 
constant ac also becomes smaller and seems to consistently decrease from above to the zero temperature result, and 
b) the critical temperature Tc increases as a fuction of the coupling constant a, with the possible saturating behavior 
in the strong coupling limit, namely the critical curve suggests the existence of the critical temperature, above which 
the chiral symmetry is always restored. 

With the data in the the critical curve we may estimate the critical coupling constant ac in the limit T ^ 0, 
ac{T — > 0), by assuming the functional form (a — adT — > 0)) oc . The fitting curve for T — > is shown in Figure 6, 
by the dashed curve, with the estimated value ac{T ^ 0) = 2.58 (l/ac(T ^ 0) = 0.388). The estimated value, 
acir ^ 0) = 2.58, is, however, significantly larger than the value ac{T = 0) = 7r/3 determined by theoretical analyses 
[1] of the DS equation for the fermion self-energy part at zero temperature, T = 0, in the ladder approximation 

in the Landau gauge with the tree level photon propagator. 

This result could be expected because of the approximation we have made use of in the present analysis. As we 
have already noticed, even if the photon propagator is HTL resummed, we adopted the improved IE approximation 
to the photon propagation. This approximation actually limits the applicability of our analysis to the rather high 
temperature region. In fact, at lower temperatures T < O.IA, the stability of the solution is getting spoiled. 
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FIG. 4: The critical coupling constant ac as a function of T. 
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FIG. 5: The critical temperature Tc as a function of a. 



We may also estimate the value of the critical temperature Tc in the strong coupling limit a oo, Tc{a oo), 
by assuming the functional form {Tc{a — > oo) — T) oc {a~^Y . The fitting curve for a — > oo is shown in Figure 6, 
by the dash-dotted curve, with the estimated value Tc{a oo) = 0.155A. The existence of the critical temperature, 
T(.{a oo), above which the chiral symmetry may always be restored, was also claimed by Kondo and Yoshida [6]. 

The critical exponents arc determined as shown in Table II. As can be seen from the table II, the estemated value 
of the exponent r] is fairly stable over the range of temperature T. The averaged value of the critical exponent r] is 



r) = 0.44 



(III.3) 



Above fact that the critical exponent t] is cleanly determined in the present analysis, strongly supports the chiral 
phase transition of QED at finite temperature under the change of the coupling constant to be of second order. As 
for the exponent ly, the estimated value is not so stable over the range of coupling constant a. The averaged value of 
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FIG. 6: The critical curve in the (T, a~'^)-planc. The dashed curve shows the fit of the critical curve for T ^ 0, and the 
dash-dotted one the fit for a — »■ oo (1/a —^ 0). The estimated values, ac{T — > 0) and Tc{l/a —^ 0) are shown with the symbol 
0, see text. 



TABLE II: Critical exponents r] and f. 



T/A 


V 


a 


V 





100 


0.450 


2.75 


0.504 





105 


0.443 


2.90 


0.562 





110 


0.452 


3.00 


0.535 





115 


0.447 


3.25 


0.614 





120 


0.436 


3.50 


0.538 





125 


0.442 


3.75 


0.469 





130 


0.445 


4.00 


0.616 





135 


0.439 


4.25 


0.780 





140 


0.425 


4.50 


0.498 





145 


0.459 


4.75 


0.780 





150 


0.446 


5.00 


0.644 



the exponent v is 



V = 0.60 . 



(III.4) 



We must improve our analysis to get stable values of v in order to understand clearly the nature of the temperature- 
dependent phase transition of QED. 



IV. CONCLUSION AND DISCUSSION 



In this paper we investigated the consequences of the improved ladder (point vertex) DS equation in QED for the 
fermion self-energy with the fully HTL resummed photon propagator, where the gauge is fixed to the Landau gauge. 
We have made use of the improved instantaneous cxchancge (IE) approximation to the photon exchange, namely, we 
have kept the exact HTL resummed transverse (magnetic) propagator, while taking the IE limit to the logitudinal 
one. This approximation can account for the essential temperature effect of the HTL resummation to the photon 
propagator. The prices we should pay due to the two additional approximations, the ladder approximation and the 
IE approximation, are the loss of the gauge invariance of the results, and the limitation of the applicability of our 
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analysis to the rather high temperature region. Several discussion and comments on these points are given at the end 
of this section. 

Before summarizing the results of the present investigation, we should stress the fact that the existence of non-trivial 
imaginary part in the fcrmion self-energy function S ^ is essentially important in the analysis of dynamical (fermion) 
mass generation at finite temperature. As we have shown in section II. 3, we cannot naively neglect the imaginary 
part without facing the inconsistent constraint equations. 

With this fact in mind, our results can be summarized as follows. 



1. The chiral phase transition is found to proceed through the second order transition. This can be seen with the 
generation of the dynamical fermion mass at the critical temperature or at the critical coupling constant without 
any discontinuity, see Figures 2 and 3. The clean determination of the critical exponents also strengthen this 
fact. 

2. The critical temperature Tc at fixed value of a is significantly lower than the previous results [6,7,8], namely, 
the restoration of chiral symmetry occurs at lower temperature than previously expected. This fact shows 
that in the previous analyses the important temperature effects were thrown away due to the inappropriate 
approximations. 

3. The critical temperature increases as a function of the coupling constant a, with the saturation behavior in 
the stronger coupling constant region, see Figures 5 and 6. This fact may indicate that the chiral symmetry is 
always restored at sufficiently high temperature, no matter how large the coupling constant becomes, agreeing 
with the result of Ref. [6] and contradicting that of Ref. [8] . 

4. The critical curve determined. Figure 6, shows that the area of the chiral broken (symmetric) phase is significantly 
smaller (larger) than previously expected [6,7,8]. This fact proves that the effect of thermal fluctuation to liberate 
the chiral symmetry is stronger than previously considered. 

5. The critical exponent 7y, that describes the phase transition at finite temperature under the change of the 
coupling constant, is determined cleanly, r\ = 0.44. This fact also supports strongly the phase transition of QED 
at finite temperature to be of second order. 

All the above results show the importance of correctly taking the thermal effects into the analysis of chiral phase 
transition. Inclusion of only the HTL resummed gauge boson (photon) propagator has significantly changed even 
the qualitative behavior of the phase transition. This fact proves a 'posteriori the effectiveness of the DS equation 
in the HTL approximation, thus also indicates the correctness of our research-strategy, namely, the importance of 
the full HTL resummed DS equation analysis of the chiral phase transition at finite temperature/density. Further 
investigation along the line of our strategy is needed to answer step by step the question how we can closely take the 
essential thermal effects into the "kernel" of the DS equation. 

As noted above, however, in the present analysis, there are two shortcomings. One is the loss of the gauge invariance 
of the result due to the ladder approximation, and the other the limitation of the applicability of our analysis to the 
rather high temperature region due to the IE approximation. The first improvement we should make is to resolve 
these problems. 

As shown in section III. 3, see Figure 1 and Table I, all the fermion wave function renormalization constants 
receive more than 10-20 percent corrections. This fact means ^ \^ and shows the explicit violation of the gauge 
invariance of the result, since Z\ = \ due to the point-vertex approximation, thus indicating the necessity of the 
gauge-parameter dependent analysis. In order to maximally respect the gauge invariance of the result in the improved 
ladder approximation, we should keep explicitly the gauge-parameter dependent contributions to the DS equations. 



Eqs. (II.8)-([I.10), and solve them with the constraint A{P) = 1, namely, Re{A{P)\ = 1 and lm{A{P)\ = 0, which 
guarantees Z2 — 1 as required, being consistent with the Ward identity Zi = Zi. This can be carried out by 
successively adjusting the gauge-parameter ^ in solving the above equations (, which may correspond to choose some 
non-linear gauge). The result of this "gauge invariant" analysis within the improved ladder approximation will be 
given in a separate paper [14]. 

The second problem of the limitation of the applicability of our analysis to the rather high temperature region due 
to the IE approximation, casts the shadow on the reliability of estimated value of aciT ^ 0) = 2.576. Actually, the 



T ^ limit of the DS equations with the IE approximation, Eqs. ([I.8)-(II.10), do not coincide with the DS equations 
in the vacuum (T = 0) theory. Thus to get more reliable estimate of the critical coupling constant in the zero- 
temperature limit, we must perform the analysis by getting rid of the IE approximation to the photon propagation, 
which is now under invetigation. 
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